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Abstract. In this article statistical bounds for certain output characteristics 
of the M /GI /1/n and GI /M/l/n loss queueing systems are derived on the ba- 
sis of large samples of an input characteristic of these systems, such as service 
time in the M/GI /l/n queueing system or interarrival time in the GI /M /1/n 
queueing system. The analysis of this article is based on application of Kol- 
mogorov's statistics for empirical probability distribution functions. 



1. Introduction 

In theoretical problems of queueing theory, the input characteristics such as in- 
terarrival and service time distributions are assumed to be known. For example, if 
we speak about an M/GI /I queueing system, we assume that the arrival process is 
Poisson with rate A, and service times are independent and identically distributed 
random variables with a given probability distribution function B{x). In practice 
we, however, have only input data characterizing arrival and departure processes, 
and our conclusion about output characteristics depend on accuracy of approxima- 
tion of the aforementioned input characteristics. 

In certain queueing systems some output characteristics can be insensitive to 
the type of probability distribution function of a service time. For example, in 
M /GI /m/Q queueing systems the stationary state probabilities are defined by the 
Erlang-Sevastyanov formulae, which are independent of the type of the probability 
distribution function of a service time. In the M/G/1 queueing system, in which 
arrival process is Poisson with rate A, and the expected service time is 6 < j, the 
expected length of a busy period is jz^i i-S- insensitive to the type of probabil- 
ity distribution B(x). Therefore, the problem of estimation of the expected busy 
period and certain other output characteristics, such as, for example, the expected 
number of served customer during a busy period, reduces to those estimations of 
the parameters A and 6 of the queueing system. 

The queueing systems, for which the output characteristics are insensitive to the 
types of probability distributions of input characteristics, are rather exceptions. In 
most cases output characteristics depend on probability distribution functions, and 
this dependence can be very complicated. So, the problem of estimating the output 
characteristics of queueing systems is generally difficult problem. 

In the present paper we demonstrate the methods of statistical analysis of certain 
output characteristics of single-server loss queueing systems such as M/GI /1/n and 
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GI lM/\/n, which are not insensitive, but expressed specificaUy via transforms of 
probability distribution function of interarrival or service times. We estimate the 
output characteristics of these systems, such as the expected busy period, expected 
numbers of lost and served customers during a busy period for the M /GI /\/n and 
the stationary loss probability for the GI /M /1/n queueing systems on the basis 
of known information about input characteristics of these systems. The concrete 
problem formulations are given later. 

Statistics of queueing systems is a distinguished area of queueing theory. The 
first publications in this area appeared long time ago (see the textbook of Ivchenko, 
Kashtanov and Kovalenko \n\ , review papers of Bhat and Rao |13j and Daley [iS] 
and the references in these sources) . In the textbook [17] various traditional meth- 
ods of statistical inference to queueing problems have been demonstrated. In [13j, a 
review of different aspects of queueing systems, including identification of models, 
parameters estimation by the maximum likelihood method and the method of mo- 
ments as well as estimates of mean value processes and auto-covariance functions, 
hypothesis testing and other topics of statistical analysis up to the publication 
date is made. In |15j . a review of various aspects, including statistical, concerning 
output or departure processes of G/G/s/N queueing systems is made. For some 
recent publications in the area of statistical analysis of queueing systems see also 
[TI] , [T^ . [IS], [21], [23] ■ Nevertheless, despite their importance, the papers on 
statistical inference of queueing systems appears much rarely compared to many 
of those that use the methods of stochastic analysis, optimization, control and as- 
ymptotic methods of Mathematical analysis. Moreover, the methods of statistical 
analysis of queueing systems known from the literature are traditional. 

The statistical analysis of the present paper is based on application of Kol- 
mogorov's statistics characterizing empirical probability distribution function of an 
interarrival or service time on the basis of large number of observations of these 
input characteristics. 

To our knowledge, Kolmogorov's statistics are never used in statistical analysis of 
output characteristics of queueing systems. In recent paper [B], Kolmogorov's met- 
ric (which is associated with one of Kolmogorov's statistics) is used for establishing 
the bounds for the loss probability in certain queueing systems with large buffers. 
In other papers 7 and [8| Kolmogorov's metric is used for establishing conditions 
for the continuity in the M/AI/\/n queueing system and, respectively, for the con- 
tinuity of non-stationary state probabilities in large closed queueing networks with 
bottlenecks. 

In the present paper we solve the following problems. Consider, for instance, the 
M /GI /1/n queueing system in which the arrival rate A and the expected service 
time h are assumed to be known, however, the probability distribution function 
B{x) of a service time is unknown. On the basis of N observations of service times 
we build an empirical probability distribution function Bcra-pix^ N). Let us denote 

6n = sup \B{x) - Bcmpix, N)\. 

x>0 

Using this statistic, we estimate the expected busy period, and other significant 
characteristics such as the expected numbers of served and lost customers during 
that busy period. We find confidence intervals (ranges) for these characteristics 
based on the confidence probability P. The similar ranges will be obtained with 
the aid of other Kolmogorov's statistics mentioned later. 
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Recall that according to Kolmogorov's theorem (see [20] or [25], p. 170) we have: 
(1.1) 



lim Pr i (5w < ^ } = K(z) 



^\EjZo.i-^ye-'' ' , forz>0, 
[0, for z < 0. 

So, for the range 

Sn < e, 

based on a chosen confidence probability P, the value e = depending on the 
large parameter N can be chosen such that the equation K{z) = P is satisfied. 
Along with 6n, the other relevant Kolmogorov's statistics are as follows: 

(1.2) 6- = sup [B{x) - Bcmp(a;, N)] , 

x>0 

and 

(1.3) S+ = sup [Bcmp(x, N) - B{x)] . 

x>0 

It is known (see e.g. [33], p. 170) that both (jl.2p and (|1.3|) have the same limiting 
distribution: 

F{z)= lim Pr|5-<^|= lim Pr { 5+ < ^} 

(1.4) r 

^ ^ _ \ l~ e"2^ ^ for z > 0, 

~ |o, for z < 0. 

For fixed N the statistics (5^ and Sj^ are dependent. However, as N increases 
to infinity, they become asymptotically independent. Explanation of this fact is 
discussed later in Section |4?2] (see Lemma 1431 and its proof). This fact is essentially 
used in our further analysis. 

For the GI/M/l/n queueing system, similar Kolmogorov's statistics are used 
for statistical analysis of stationary loss probabilities, where the probability dis- 
tribution function of interarrival time is denoted by A{x), empirical probability 
distribution of interarrival time based on TV observations of interarrival times is de- 
noted by Aerap{x,N), and statistics Sn, SJ; and S'j^ are, correspondingly, as follows: 

Sn = sup \A{x) - Acnip{x, N) I , 

x>0 

(5^ = sup [A{x) - ^cmp(a;, N)] , 

2;>0 

and 

Sjf = sup [^cmp(a::, N) - A{x)] . 

x>0 

The aim of the present paper is twofold. Along with obtaining the ranges for 
aforementioned characteristics of M/GI/l/n and GI/M/l/n queueing systems for 
a given probability P (which is the first aim of the paper) we also answer to the 
question: Which of these Kolmogorov's statistics is better? That is, under which 
of these Kolmogorov's statistics the difference between an upper and lower bounds 
based on the chosen probability P is smaller? 

The rest of the paper is organized as follows. In Section 2, we recall the known 
representations for the characteristics studied in the paper for M/GI/l/n and 
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GI /M/\/n qucucing systems. In Section 3, we prove the main lemmas, which are 
then used to find the estimators for required characteristics of the queueing systems. 
The estimators themselves are derived in Section 4. In Section 5 we compare our 
statistical results and address the question formulated above. In Section 6 numerical 
results are presented. In Section 7 we conclude the paper. 



2. The recurrence relation for main characteristics of the 
M/GI/l/n AND GI/M/l/n queueing systems 

In this section, we recall the known results for the main characteristics of the 
M/GI/l/n and Gl /M /\/n queueing systems during their busy periods. For a 
more detailed information see [T], [5] or [TU]. 

2.1. The M/GI/l/n queueing system. Consider the M/GI/l/n queueing sys- 
tem, in which the arrival flow is Poisson with parameter A, and the probability 
distribution function of the service time is B{x) having the expectation b. Param- 
eter n denotes the number of waiting places, i.e. the capacity for the customer 
in service is not taken into account. Let r„ denote the length of a busy period 
of this system, and let Un and L„ denote the number of served and, respectively, 
lost customers during that busy period. The recurrence relation for ET„ has been 
originally obtained by Tomko |26j : 



(2.1) ET„ = ^Er„_,+i / e-^-^di3(x). 



where ETq = 6, and Ti denotes the length of a busy period in the M/GI /1/i queue- 
ing systems having the same arrival rate A and the same probability distribution of 
the service time as the original queueing system M/GI/l/n. The expected number 
of served customers during the same busy period satisfies the recurrence relation 
similar to (12.11): 



(2.2) E;.„ = ^Em„_,+i / e-^-^-^dB{x), 



where i^i denotes the number of served customers during a busy period Ti (e.g. see 
[2]). The main difference between recurrence relations (|2.ip and (|2.2p is that (|2.2p 
starts from Ei/q — 1, while p.ip starts from ETq — b. 

For the expected number of losses during a busy period we correspondingly have 
the following recurrence relation (see e.g. [I], [2]: |S] or [TU] for more details), which 
is similar to the previous two recurrence relations given by (|2.ip and (|2.2p : 

(2.3) EL„ - 1 = V(EL„_,+i - 1) r c-^-^dB{x), 

'H. Jo 



j=0 



where ELq — ^b, and Li denotes the number of losses during a busy period Ti. 
All (|2.ip . (|2.2I) and (12.31) are the following convolution type recurrence relation: 



(2.4) Qn = 



i=Q 



where = /o°° e Ax (Ax) ^j^j^^^y 
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2.2. The GI /M /1/n queueing system. Consider the GI/AZ/l/n— 1 queueing 
system, where an interarrival time has the probabihty distribution function A[x), 
and the parameter of an exponentiaUy distributed service time is /i. (The number 
of waiting places n — 1 excludes the place for a customer in service.) Let 7r„ denote 
the stationary loss probability. It is shown in 3 that 7r„ satisfied the recurrence 
relation: 

(2.5) re--i^dA(.), 

where ttq = 1. Recurrence relation (12. 5p is similar to that of (|2.ip . (|2.2p and 
(|2.3p . and has general form (|2.4I) . The only difference is that in representation 

(^3]) r, = /Q°°e-'"^i^dA(a;), while in representations (EI]), dO]) and (ES]) r, = 

re-^^^dB(x).^' 

Another representation for 7r„ has been obtained by Miyazawa [35]. 



3. Main lemmas 

In this section we prove main lemmas, which help us to obtain then the desired 
estimators for characteristics of queueing systems studied in this paper. 

Let Fi (x) and f 2 {x) denote arbitrary probability distribution functions of posi- 
tive random variables. For a positive parameter a, let us denote 

Jo 



and 



r^(F2)= / e 
/o 



7.1 



Lemma 3.1. Assume that snp^y^^ ~ P2{x)\ < e. Then, 

(3.1) |ro(Fi)-ro(F2)| <e, 
and for all i — 1,2, 

(3.2) |r,(Fi)-r,(^^2)| <2e. 
Proof. By partial integration we have 

POO POO 

(3.3) / e-"^dFfc(x) = a / e-"^i^fc(a;)dx, fc = 1, 2. 
Jo Jo 

Therefore, 



/ e-"Mi^i(x)- / e-"^dF2(a;) =a / 
Jo Jo Jo 



e-"=^i^i(x)da; 



e-"^F2(.T)da: 



<a [ e-"^da: (sup | Fi {x) - F2 {x) 

Jo \x>0 



<e by the assumption 



< e. 



Inequality p.ip is proved. 
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Let us now prove inequalities p.2|) . For i = 1,2,... and k = 1,2 by partial 
integration we have: 



(3.4) 



-ax 



(ax) 



« 

oo 

a [ I e 



[ax 



-Fk{x)dx 



Therefore, 



(3.5) 



7,' 







■^^dF.ix) 



^(^Fi(x)da;- / e 
«! Jo 



,{ax) 



-^F2{x)dx 



< 



^— ^Fi(.)da: + y^ e ^— ^F.(:.)dx 
/ e-"="(aa;)Mx(sup|i^i(x)-F2(x)| ) 

" v ' 

<e by the assumption 



(*-l)Uo 



(ax)*"Mx sup \Fi{x) -F2ix)\] . 



x>0 



by the assumption 



Notice, that 



and 



/ e-"^(ax)Ma; = - / e'yy'dy = -r( 
Jo a Jo " 



i + 1), 



/ e-"^(ax)'-M2; = - / e'yf-^dy = -r(i), 
Jo "Jo " 

where r(a;) is Euler's Gamma function. Taking into account that T{i + 1) = «!, 
from the last inequality of (|3.5p we arrived at the desired result \ri{Fi) — ri{F2)\ < 
2e for i = 1,2,.... Inequalities p.2p are proved, and the proof of the lemma is 
completed. □ 



Lemma 3.2. Assume that sup^^^Q [Fi{x) — -^2(2^)] < ei and sup^-^g [^2(2;) 
Fi(x)] < £2- Then, 



(3.6) 



roiFi) - ro{F2) < ei. 



and for all i — 1,2, 



(3.7) 



n{Fi)-n{F2) < €1 + 62. 
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Proof. The proof of this lemma is similar to that of Lemma 13.11 From p.3p we 
have: 

POD fOO fOD POD 

/ e-"MFi(a;)- / e-""dF2(x) = a / e-""Fi(x)da; - / e-""F2(a;)dx 
Jo Jo Jo Jo 

<a f e-°"'dx (sup [Fi{x) - i^2(x)] 

Jo \x>0 


<£i (assumption) 

< ei. 

Inequality p.6p is proved. Next, from (|3.4p we have: 

e-"^-^d^^i(a;) - H e-"^^^dF2{x) 
Jo ^ 

a /"e-""-^i^i(a;)dx- e-"^^^F2(a;)dx 



< 



Jo 

/•OO / 

- / e-"%axydx sup [^1(2;) - F2{x 

»' Jo V:r>0 

=^ <ei (assumption) 

+ / e-""(ax)'-Ma: fsup [F2{x) - Fi{x)]] 

v~ » ' 

= <e2 (assumption) 

< ei + £2. 

Desired inequality p.7p follows. The proof is complete. □ 
4. Explicit recursion and estimates for characteristics of queueing 

SYSTEMS 

The recurrence relation in form (|2.4p is not convenient for estimating the afore- 
mentioned characteristics of the queueing systems. In the left-hand side of (12. 4p 
Qn is, while in the right-hand side of (12.41) the linear combination of Qi, Q2, - ■ ■ , 
Qn+i is presented. In order to derive the appropriate estimates, we first should 
rewrite (|2.4p in the form of explicit recursion. We have: 

(4.1) Qi = -Qo, 

ro 

(4.2) Q2 = ^—^Qi: 

ro 

(4.3) Qn = — [{I ~ ri)Qn~i - r2Qn-2 - ■ ■ ■ - rn-iQi] , n >3. 

ro 

Explicit recurrence relations of (|4.ip - (|4.3l) can be now used to establish necessary es- 
timates for characteristics of the loss queueing systems M/GI/l/n and GI/M/l/n. 
To build the recursion for ET„, we start from the M/GI/1/0 queueing system. The 
busy period of this system contains only a single service time. In the sequel we as- 
sume that the parameters A and b of the queueing systems are known. As the 
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expected service time b is known, we set To = b, that is, the expected busy period 
of the M/GI/l/0 queueing system is reckoned to be estimated exactly. In the 
foUowing we use the notation: 

n{B)^ re-^-^di3(x), 
Jo 

and 

r,(Somp)= /"e-^-^dBemp(x,iV), 

where Bcmpi^, N) is an empirical probability distribution function based on N ob- 
servations. In the following the estimator for Q„ will be denoted Qn and specifically 
that for the expected busy period will be denoted Tn- (Similarly, for GI/M/l/n 
queueing systems, the notation ri{A) and ri(Aonip) can be used where the parameter 
A should be replaced with the parameter ^.) 

4.1. Estimators based on the statistic Sj^. Assume that the inequality 
(4.4) Sn < e 

holds with probability P. Since N is assumed to be large enough, this probability 
can be chosen from limit relation p.ip . 

Then, according to Lemma l3.1i with the probability not smaller than P we have 

|ro(-B) - ro(-Boinp)| < e, 

and for z = 1, 2, . . . , n 

\ri{B) ~ ri{B^n-ip)\ < 2e. 

For further simplifications, we will write = f j(i?cmp), « = 0, 1, 2, . . ., omitting the 
argument i?cmp hi the notation. 

On the basis (|4.ip - (|4.3p and Lemma [01 we have the following. 

Theorem 4.1. The point estimator Qn for a required characteristic of a queueing 
system is recurrently defined as 



(4.5) 


Qi — — Qq, 






(4.6) 


Q2 — Qi, 

ro 






(4.7) 


ro 


(1 - ri)Qn-l - r2Qn-2 - ■ ■ 


• - rn^iQi 


, n > 3 



Then, the interval estimators with a confidence probability non-smaller than P are 
recurrently defined as 

(4.8) qT"- = -^Qo <Qi< ^— go = gi^p"', 

J-o + e ro-e 

(4-9) g^- - \7;'f go <Q2< \7iy go = gr''^ 
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■flower 



(1 - n - 2e)Q^— 



(4.10) 



-{r, + 2e)g™ - ... - (r„_i + 2e)Q^PP" 

< Qn < f(l -n+ 2e)QZT 

ro - e L 

-(r2 - 2e)gJ°-^'- - ... - (r„_i - 2e)Q'r-" 
= Q™, n > 3. 

T/ie uaZfie e is determined from the equation: K (^e\/jfj = P, ^^(2;) is Kolmogorov's 

Junction defined in (II. ip . and m t/ie notation of ()4.8p - (l4.10p Q'"™'' and Q"^^" are 
itsec! /or lower and, respectively, upper bound of the interval estimators. 

Proof. Indeed for any < e < tq we apparently have 



(4.11) 



lower 



1 



7-0 7-0 - e 



_ ^uppor 



where -^Qo in the middle of inequality (|4.1ip is Qi. This means that for e > given 
such that |ro(S) — ro(-Bcmp| < e occurs with confidence probability not smaller than 
P, we have 

qT" = ^— Qo < Oi < ^— go = Qi"^p". 

ro + e ro - e 

Relation (|4.9p follows similarly. Indeed, for any < e < ro we apparently have 



(4.12) 



l-ri-2e 1-^1^ l-ri+2e 



[^0 



where ^-^Qo in the middle of inequality (I4.12p is Q2- That is, for any e > given 
such that |ro(i?) — ro(i?cmp| < e occurs with the confidence probability not smaller 
than P, we have 

l-n~2e l-r,+2e 
[ro + 6]2 ^° < O2 < Qo. 
Then (|4.10p is proved by the similar arguments with the aid of induction. Namely, 
for n > 3 and any < e < ro we have the inequality 

^ [(1 - n - 2e)te^ - (r2 + 2e)g™ - ... - (r„_i + 2e)Q\^'''A < Q„ 



ro + e 



< 



1 



ro - e 



(1 - n + 2e)Q:r_r - (r2 - 2e)Q'°-l 



(r„_i - 2e)g'/ 



lower 



which means that for e > given such that |ro(P) — ro(Pomp| < e occurs with the 
confidence probability not smaller than P, we have 
1 



(1 - n - 2e)Q;,°-- - (r2 + 2e)QZT (rn-i + 2e)Qf 



< Qn 



< 



1 



;i - n + 2e)QZT - ir2 - 2e)Qjr^'- 



ro - e 

Theorem 14. 11 is proved. 



(r„_i - 2e)g'°' 



□ 



Using relations (I4.5p - (|4.10p one can build estimators and the bounds for these 
estimators for all of the aforementioned characteristics of the queueing systems 
M/GI/l/n and GI/M/l/n. 
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Remark 4.2. For relations (|4.8p - (|4.10p the following conventions are made. If ro < e, 
then the right-hand sides of (|4.8p - (|4.10|) are set to be equal to infinity. If 1 — ri — 2e 
is negative, then the left-hand side of (|4.9|) is set to be equal to zero. In addition, 
if one or other term 1 — — 2e, i = 1, 2, . . . , n of the left-hand side of (I4.10p is 
negative, that term is set to zero. The value e should be chosen small enough, i.e. 
appropriate sample size N should be chosen large enough, such that all of the above 
inconsistences can occur with negligibly small probability. 

4.2. Estimators based on the statistic or 5^. Since both of these statistics 
are similar, we only consider the first one, 6~^. Assume that N is large enough and 
6^ < e occurs with the probability not smaller than P. 
According to Lemma 3.2 with the same probability P 

ro{B) - ro(Bcmp) < e, 

and 

ro{Bcmp) - ro{B) < e. 

From same Lemma [3T2] according to p.7p we have r-,;(i?) — ri(_Bcmp) < + ^o^' 
i — 1,2, ... ,n. Therefore, given that Pr{(5^ + 5^ < 7} = P, with the confidence 
probability not smaller than P, we obtain the inequalities: 

ri{B) - ri(Bcmp) < 7i 

and 

n{B^mp) - n{B) < 7. 

Lemma 4.3. As N increases to infinity, the statistics y/NS^ and y/NS'j^ become 
asymptotically independent. 

Proof. This fact can be proved with the aid of construction used in [22, p. 171-172 
as follows. Assume that B{x) is a continuous probability distribution function, 
and Bcmpix, N) is the empirical probability distribution function based on the 
sample (^1, £,2,- ■ ■ ,£,n), that is, Bempix, N) is defined as the fraction of random 
variables that less than or equal to x. Let us write the sample (^1, £,2,- ■ ■ ,£n) 
in ascending order as {£1, £,2t--j£n)- Consider the deviations S'^{n) — B{£^) — 
B^^p{C,N) and 6+{n) = Pomp(C*,iV) - B{C), n = 1,2,...,N. Apparently, 
the random variables S'^{n) and S^{n) all are continuous (n = 1, 2, . . . , N), and all 
their joint distributions are independent of the choice of the probability distribution 
B{x). That is, without loss of generality one can assume that B{x) is the uniform 
distribution in [0,1]. Then, Sj^{n) = - and Sj^{n) = 77 - C- Let k„ be the 
number of random variables amongst ^1, £2,- ■ ■ , £n falling into the interval (^^^^, 
and Afn = ki + K2 + . . . + k„, n — 1,2, . . . , N . Then the number of points where 
S'^{n) are positive is X]^=i I{-^n — 3t.nd, respectively, the number of points 
where are positive is ^n=i^i-^n < n}. For N fixed, the random variables 

— maxi<„<Ar i5^(n) and S'j^ — maxi^nKN S^in) are generally dependent. As 
N — )■ cxD, the random variables and 5^ both vanish. The aim is to prove 
that the normalized random variables \/N6^ and ^/NS^ become asymptotically 
independent as — > 00. 

For this purpose we will use the known fact that the random variable ^* has the 
Beta distribution with parameters n and A^ — n + 1 (see e.g. David and Nagaraja 
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[l4] . p. 14). In other words, this means that the random variable ^* has the same 
distribution as 

7?1 + 7?2 + ■ ■ ■ + 
111+112 + ■■■+ VN+l ' 

where rji, 772,- ■• , Vn,- ■ - is the sequence of independent exponentially distributed 
random variables with parameter 1 (see Karlin [TH], p. 242-244). According to the 
strong law of large numbers, 

Prj lim — = ll = l. 

[W^oo 771 + 772 + • ■ ■ + W+l J 

Therefore, as iV — )■ 00, the random variable NSj^{n) = N^*^ — n can be asymptoti- 
cally represented as 

(7yi-l)-K(772 -l) + ... + (77„-l). 

From this asymptotic representation, apparently that limAr_j.oo E(A^^* — n) = for 
all n = 1, 2, . . . , A^, and the expected number of positive and negative values of 
A^* — n for n = 1,2, ... ,N is approximately the same. That is the numbers of 
points where S~l^{n) and Sjf{n) are positive are asymptotically equal. 

For large A^, let j denote an index such that NS^{j) is negative while A^i5^(j-|-1) 
is positive. Then there is a random number v of consequent positive values of 
NS^{j + i), i = 1,2, . . . ,v. This random number is called length of positive period 
or simply positive period. Notice that, 

(4.13) Pr <^ ^(r;, - 1) < x| ^(7;, - 1) < 0, ^(77, - 1) > = 1 - e'^ 

lz=l 1=1 i=l ) 

Indeed, let Tj = j — X]i=i Vi + ^- Then, the conditional probability of the left-hand 
side of (|4.13l) can be written 

/>oo 

Pr{7/j+i <x + Tj\r/j+i > Tj > 1} = / Pr{7/j+i <x + y\r]j+i > y}dPr{r, < y} 

Jq 

1*00 

= (1 - e--) / dPr{T, < y) 
Jo 

= 1-e-^. 

(j4.13p follows. Thus, the limiting distribution of NSjf{j + 1) as A^ — 00 is ex- 
ponential, and NS~l^{j + 1) is asymptotically independent of the past values of 
NS^{k), k = 1,2, ... ,j. The following random variables N6^{j + i), i ~ 2,3, . . . ,v 
are also asymptotically independent of the aforementioned past values NS^{k), 
k = 1,2, ... ,j. So, the positive period and associated random variables NS^{j + i), 
i = 2, 3, . . . , u all are asymptotically independent of the past. 

Consequently, let ji, j2,- ■ .he the sequence of values such that N6^{jk), k = 
1,2,... are negative while NS'^{jk + 1) are positive. Then, according to the above 
statement, the positive periods are asymptotically mutually independent and iden- 
tically distributed random variables. 

The similar result can be obtained for negative periods which are constructed 
similarly to those positive periods. Indeed, for A^ large, let now j denote an index 
such that N6^{j) is positive while N6^{j + 1) is negative. Then there is a random 
number t of consequent negative values of NSj^{j + i), i = 1,2, . . . , l. This random 
number is called length of negative period or simply negative period. Notice that 
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the event {NS^{j) > and NS^{j + 1) < 0} means that during the time interval 
[j, j + 1) there are at least two events of Poisson process with rate 1, and the second 
event can be considered as the start of a negative period. As in the case of positive 
period, the distribution of the length of a negative period does not depend on the 
past history as well. 

Thus, as TV — > oo, positive and negative periods, alternatively changing one an- 
other, are asymptotically mutually independent and identically distributed random 
variables. 

As iV — >■ oo, the sets of values where N6^{n) and NS]^{n) are positive, are 
asymptotically independent, because the positive and negative periods are asymp- 
totically independent. Note, that the set of indexes where the periods are positive 
is {n < iV : Afn < n} and that set where periods are negative is {n < N : Mn > n}. 
Furthermore, 



Pr I lim y I{AA„ < „} = i I = 1 



On the other hand, as — cx), the limiting distributions of \/lV6^ and \/a'5^ 
do exist, and they are defined by (|1.4p . 
Therefore, 

lim Pr (\/iV(5+ < X, ViV(5jv < y] 

= lim Pr i max N6tj(n) < x, —= max NS7r(n) < y 

N^oc [ y/N {n<JV:A/'„<n} ^/N {n<JV:A/'„>n} 

= lim Pr 1 —= max N6tr(n) < x 



X lim Pr i —;= max NSmM < y 

N^oo 1 ^/]V {n<N:Ar„>n} 



lim Pr (y7V5+ < x] lim Pt\VN6]^ < y\ 



and the proof of the statement of the lemma is completed. □ 
From Lemma 14.31 and (|1.4p we have: 

(4.14) ^lim^Pr|<5^ + <5+<-^|=i^*i^(z) 



where the probability distribution function F(z) is defined in (11.41) . and the asterisk 
denotes convolution of the probability distribution function F{z) with itself. 

For the convolution of the probability distribution F{z) = 1 — e~^^ with itself 
we have: 

(4.15) / (l-e-^(^-^)')4a;e-2^'dx==l-e-2^'-V^ze-^'[2$(y2z)-l], 
Jo 

where ^(z) = J^^e~^ dy. So, the right-hand side of (j4.14p is explicitly de- 
termined. 

Now, we have the following theorem. 

Theorem 4.4. The point estimator Q„ for a required characteristic of a queueing 
system is recurrently defined by (|4.5p - (|4.7p . Then, the interval estimators with a 
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confidence probability non-smaller than P are recurrently defined as 



(4.16) 
(4.17) 



1 



flower 



1 



■To 



-Qo <Qi< 



1 



:Qo — Qi^ 



lower 



[^0 + e] 
1 



ro - e 



1 - ri + 7 



^ upper 



(1 - n - 7)QLTr 



(4.18) 



-(r2+7)Qrr 
1 



< Qn < 



ro 



flower 



..-(r„-i+7)Qr'=^ 
■ - (rn-i - 7)Ql°""' 



= or^"', " > 3, 

where the value e is determined from the equation: 1 — exp (— 2A^e^) — P, and the 
value 7 is determined from the equation: 

1 - exp {-2N-f^) - V^7 exp (-iV7^) 2$ (y2N-f^ - 1 = P, 
where $(x) — e~^da;. 



Proof The proof of (|4.16p - (|4.18p is similar to that proof of 



5. Which of the statistics is better? 



(14.101) in Theorem 
□ 



In this section we address the question which of the statistics 6n, or Sjf is 
better, or under which of them the difference between upper and lower bound is 
smaher? The statistics 5^ and Sjj are symmetric and they give the same bounds. 
Therefore the question of comparison should be addressed to the statistic Sn and 
one of the statistics or (5^, say (5^. 

Clearly, that 



(5.1) 



Sn := sup \B{x) - Bcn-,pix,N)\ 

x>0 

>sup [B{x)- B^n,p{x,N)] 

x>0 



■.= s 



N- 



Therefore, for the same probability P the z-value of the equation P — K{z) is not 
smaller that the z-value of the equation P = 1 — e~^^ for any given probability 
P. As a result, ^{5^) ^ £(<5jv)' '^tiere ^{Sn) is the value e obtained for the statistic 
(5 AT, and e{5jf) is the value e obtained for the statistic 6jj both for the same value 
of probability P. 

On the other hand, for any probability distribution function G{x) of a positive 
random variable we have as follows. Let Fi{x) — G (|) and ^2(2^) = (G * G){x), 
where the asterisk denotes convolution. Apparently, 



(5.2) 



xAFi{x) 



xdF2{x), 
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and 
(5.3) 



Jo Jo 



(the integrals in (I5.2|) and (I5.3|) are assumed to converge). Then (see e.g. 
F2{x) is said to be smaller than Fi{x) in the convex sense. 

Then, there exists a point xq where the probability distribution functions Fi{x) 
and F2{x) cut one another (see Karlin and Novikoff [19] or Stoyan [2^, p. 13). 

In the given case F{x) = 1 — e^^^ . Therefore, Fi{x) = 1 — e~5^ , while 
^2(0:) = 1 — e~^^ — -y/TTxe^^ [2$(-\/2a;) — 1] ■ Apparently, there exists a point xo such 
that Fi{z) < F2(z) for all z > xq. Numerical calculations show that xq ~ 1.385, 
which corresponds to the level of probability Fi(1.385) « ^2(1. 385) ~ 0.6166. 

Let zi and Z2 be two z-points associated with the solutions of the corresponding 
equations Fi(z) = P and ^2(2) = P, where P is some given level of probability. 
Clearly, that if both zi > xq and Z2 > xq, then, since Fi{z) < F2{z) in this set, we 
also have zi > Z2- 

This enables us to conclude, that correspondingly to these z-points of the afore- 
mentioned probability distribution functions we obtain 2e{5jj) > j{Sjf). The fol- 
lowing notation is used here. As before, by £{S^) we mean the value e associated 
with the solution 1 — e^^^*^' = P. By 7((5jy) we mean the value 7 associated with 
the solution 



1 - exp (-2A^7^) - VnN'fexp (~Nj 



2\ 



= P. 



The aforementioned relation 2e{5^) > j{6^) follows from (15. 2p and (|5.3p as follows. 
For a given z value associated with large level of probability P we have Fi(z) < 
F2{z). Hence, as it was mentioned before, from Fi(zi) = P and -^2(^2) = we 
obtain 2:2 ^ ^i- That is, from (15. ip we have e((5^) < e{Sj^) and, for large enough 
z-values taking into account that F{z) ~ Fi{2z), we have "f{6]^) < 2e(d^), and 
since e((5^) < e{SN) we finally have 7(i5^) < 2e((5Ar). (Recall that by e{5N) we 
mean the value of e associated with the solution K{f\fN) = P.) 

Thus, the statistic 5jj becomes better than the statistic 5n for large enough z- 
values, i.e. for large enough values of probability P. This conclusion is supported 
by the numerical examples in the next section. 



6. Numerical examples 

Numerical examples are provided for N — 10, 000. For simplicity we take B{x) = 
1— e"''^, and A = ^ = 1. The value n is taken 4. In the numerical examples below we 
build estimators for the expected busy periods of the M /G//1/ n queueing systems. 

6.1. The statistic 5m- From the equation K{z) — P — 0.95 we obtain z w 1.3581. 
Therefore, 

1.3581 



= 0.013581. 



V10,000 

With the aid of simulation, we build empirical probability distribution, and 
on the basis of that empirical probability distribution we obtained ri(_Bcmp), i — 
0, 1,2,3,4, given in Table 1. 
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Expected 


The value obtained 


Statistics 


(theoretical) 


by simulating 




value 


empirical distr. 




r^{B) 


(-^cmp) 


ro 


.5 


.5031 


n 


.25 


.2488 




.125 


.1234 


r-3 


.0625 


.0615 


r4 


.03125 


.0308 



Table 1. The table for statistics ri{Bcmp) obtained by simulating 
empirical probability distribution function 





Expected 


Estimated 


Lower 


Upper 


Statistics 


(theoretical) 


value 


bound 


bound 




value 










ET, 




T^lowcr 
i 


flipper 

i 


To 


1 


1 


1 


1 


Ti 


2 


1.987589 


1.935434 


2.042817 


T2 


3 


2.967558 


2.71285 


3.248177 


T3 


4 


3.943322 


3.206328 


4.783615 


Ti 


5 


4.916821 


3.317455 


7.057548 



Table 2. Estimators for expected busy periods in AI/GI/l/n 
queueing systems based on the statistics 



Then, on the basis of statistics ri{Bcmp), i = 0il:2,3,4 the lower and upper 
bounds for estimators Ti, i = 1,2,3,4 and these estimators themselves are shown 
in Table 2. 

6.2. The statistic (5^ or S^. From the equation 1 — e~^^^ = P = 0.95 we obtain 
z — 1.224. Therefore e — 0.01224. From the other equation 

1 - e-^^" - V^ze^^' [2$(\/2z) - l] = P = 0.95 

(see (|4.15p ). we have z — 2.08. Therefore 7 — 0.0208. In this case the lower and 
upper bounds for estimators Ti, i = 1,2,3,4 are shown in Table 3. (For convenience 
the estimators Ti, i = 1,2,3,4 themselves are duplicated from Table 2.) 

7. Discussion and the future work 

We provided statistical analysis of certain output characteristics of M/GI/l/n 
and GI/M/l/n queueing systems by using the known statistics associated with 
empirical distribution. It follows from our results and numerical analysis that in 
certain cases the use of the statistic 6jj or 5^ is more profitable than that of 
Sjf. Namely, in these certain cases the differences between the upper and lower 
bounds for the output characteristics are smaller when we use the statistic (5^ or 
(5]^ compared to these differences between the bounds when we use the statistic 6n- 

It should be noted, that presented methods can be used for other queueing 
models, in which characteristics can be expressed via convolution type recurrence 
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Expected 


Estimated 


Lower 


Upper 


Statistics 


(theoretical) 


value 


bound 


bound 




value 










ET, 


T^ 


T^lowcr 
i 


^uppcr 
i 


To 


1 


1 


1 


1 


Ti 


2 


1.987589 


1.940466 


2.037241 


T2 


3 


2.967558 


2.750256 


3.204070 


T3 


4 


3.943322 


3.327933 


4.633603 


Ti 


5 


4.916821 


3.616202 


6.673106 



Table 3. Estimators for expected busy periods in M/GI/l/n 



queueing systems based on the statistics Sj^ or 



relation (|2.4p . An example is the extended M/GI/l/n loss queueing system de- 
scribing models of telecommunication systems considered in j4| . Another example is 
the state-depending queueing system describing models of dam/inventory systems 
considered in [S]. More complicated models of loss queueing system with batch 
arrivals/departures of customers are of special interest as well and can be consid- 
ered as a subject of future work. It is also interesting to study the cases where 
the parameters A and/or h are unknown. This case is more realistic for practical 
applications and leads to new challenging aspects of this theory. 
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Appendix: Calculation of the integral in (|4.15p 

We have: 

/ (l-e-2(^-^)')4a;e-2^'da;= / 4a;e~2^'dx- ( 4xe-2^'e-2(^-^)'da; 
Jo Jo Jo 

= h + h. 
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Clearly, that Ii = F{z) = 1 — e . For I2 we have the following: 

l2 = - f 4x6-2^'+'*^^-^^' da; = -e"^' / 4xe-(^-2^)'dx 
Jo Jo 

= -e"^' / [(2a; - z) + z]e-<^''-^''^^ d[2x - z] 
Jo 

-z J — z 



=0 



Hence 



1 /■+\/2z 

= -x/¥0e-^^ / e-* /2df 

V27r J-ySz 

= -v/^r^e"^' 2$ - 1 . 

^ (1 - e-2(^-^)')4a;e-2^'da; = 1 - e'^^' - s/^ze''^ [2$ (\/2^) - 
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